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Absence of self-averaging and of homogeneity in the large scale 
galaxy distribution 
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Abstract. - The properties of the galaxy distribution at large scales are usually studied using 
statistics which are assumed to be self-averaging inside a given sample. We present a new analysis 
able to quantitatively map galaxy large scale structures while testing for the stability of average 
statistical quantities in different sample regions. We find that the newest samples of the Sloan 
Digital Sky Survey provide unambiguous evidence that galaxy structures correspond to large 
amplitude density fluctuations at all scales limited only by sample sizes. The two-point correlations 
properties are self-averaging up to approximately 30 Mpc/h and are characterized by a fractal 
dimension D — 2.1 ± 0.1. Then at all larger scales probed density fluctuations are too large in 
amplitude and too extended in space to be self-averaging inside the considered volumes. These 
inhomogeneities are compatible with a continuation of fractal correlations but incompatible with: 
(i) a homogeneity scale smaller than 100 Mpc/h, (ii) predictions of standard theoretical models, 
(iii) mock galaxy catalogs generated from cosmological Nbody simuations. 



Introduction. — Understanding the large scale struc- 
ture of the universe as mapped by galaxy distribution rcp- 
C^fesents one of the cornerstones of modern cosmology. It 
provides the basic test for theories of structure formation 
in the universe. A primary question in the statistical anal- 
ysis of three-dimensional galaxy catalogs (where, in addi- 
tion to the angular coordinates, the redshift is measured 
and through Hubble's law [1] the distance of each object) 
concerns the determination of a scale where the distribu- 
tion becomes homogeneous. Such a scale A can be defined 
to be the one beyond which counts of galaxies in three di- 
mensional spherical volumes of radius r grow as r 3 [2]. 

A decade ago, by measuring the conditional density, i.e., 
the local galaxy density seen by a galaxy in a spherical 
volume of radius r around itself [2,3], some of us found that 
galaxy correlations are power-law with an exponent 7 « 1 
up to the sample sizes, i.e., ~ 30 Mpc/h 1 , corresponding 
to a fractal dimension D = 3 — 7 ~ 2 [4-6] . These results 
were in contrast with the analysis of the same samples by, 



1 We use H = W0h km/scc/Mpc, with 0.4 < h < 0.7, for the 
Hubble's constant 



e.g., [7-9], who found A « 10 Mpc/h and 7 = 1.8. The 
reason for these differences lies in the a-priori assumption 
of homogeneity, inside a given sample, of the standard 
statistical analysis [2-4]. 

At larger scales, with weaker statistical significance, 
there was an evidence compatible with the fact that power- 
law correlations in the conditional density extend up to 
r ~ 100 Mpc/h or more [4,10]. These results generated 
a debate in the field [11-13] because even though galaxy 
structures were found in many different catalogs to ex- 
tend to scales of the order of hundreds of megaparsecs, 
the characteristic length scale Ao statistically describing 
their correlations was determined to be a few megaparsecs 
[7-9]. While for some this was a paradox [3-5], for oth- 
ers [7-9,11,12] the explanation was that large scale struc- 
tures have small amplitude relative to the average density. 
However this interpretation is problematic as in the range 
of scales where the conditional density shows power-law 
correlations the sample density is not well defined while 
density fluctuations have large amplitude [2-4] . The deter- 
mination of the crossover scale Ao, where the conditional 
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density from a power-law turns to a constant, allowing a 
meaningful determination of the average density, has been 
thus an important task of galaxy correlations studies in the 
last decade [14]. 

Two new galaxy catalogs, the Sloan Digital Sky Survey 
(SDSS — [15]) and the Two degree Field Galaxy Red- 
shift Survey [16], have recently provided great advances 
in the mapping of the local universe both for the num- 
ber of objects measured in continuously growing volumes 
and for the determination of several parameters for each 
of them. Several studies [17-19,22] of different samples of 
these surveys confirmed the small scale correlations mea- 
sured by [4,6]. In addition it has been claimed that a slow 
crossover toward homogeneity occurs [17] with the average 
conditional density in spheres at ~ 20 Mpc/h having twice 
the amplitude of the asymptotic density reached at r > 70 
Mpc/h [18]. It was however noticed that galaxy structures 
could bias the determination of correlations in these sam- 
ples introducing uncontrolled systematic effects [19,22,23]. 
Recently in the Two degree Field Galaxy Redshift Survey 
it has been found that [20,21] galaxy distribution charac- 
terized by large 1 amplitude fluctuations with a large spa- 
tial extension, whose size is only limited by the sample's 
boundaries. In addition at scales r < 40 Mpc/h, it has 
been observed a well defined and statistically stable power- 
law behavior of the average number of galaxies in spheres 
in agreement with previous determinations. 

A different and complementary method to characterize 
structures is provided by galaxy counts as a function of the 
radial distance from us or of the apparent luminosity [24] . 
These show large fluctuations around the average behavior 
both in redshift [25] and angular surveys [26,27]. There 
have been controversies as to whether these are due to 
real clustering or to incompleteness of the catalogs [11,24]. 
Recent results support the conclusion that the local galaxy 
distribution is characterized by large scale structures with 
significant correlations on scales r > 50 Mpc/h [28,29]. 

In this paper we use a new method which is able to 
establish, in a given sample, the link between the small 
scales r < 30 Mpc/h correlations and the large scales 
r > 30 Mpc/h fluctuations in galaxy counts and which 
clarifies how the latter influence the determination of the 
former. Using it we can test whether sample means, vari- 
ances and correlations are well defined, i.e., whether they 
are statistically stable in different sub- volumes of the given 
sample. By applying this method to the data of the SDSS 
project [30] we detect large density fluctuations of spatial 
extension limited by the samples' sizes. We show that 
these introduce systematic biases in the determination of 
large scale correlations. 

The Data. - The SDSS [15] is currently the largest 
spectroscopic survey of extragalactic. Here we consider 
the main galaxy (MG) sample (DR6) [30] containing red- 
shifts for about 800,000 galaxies which covers an area of 
7425 square degrees on the celestial sphere. To query the 
DR6 database we constrain the flags indicating the type 



of object so that we select only the galaxies from the MG 
sample. We then consider galaxies in the redshift range 
10 -4 < z < 0.3. The redshift confidence parameter is 
constrained to be z con f > 0.35 with flags indicating no 
significant redshift determination errors. In addition we 
apply the filtering condition m r < 17.77, using Petrosian 
apparent magnitudes in the r filter which are corrected for 
galactic absorption, and thus taking into account the tar- 
get magnitude limit for the MG sample in the DR6 [31]. In 
this way we have selected 479,417 objects. We considered 
also more stringent limits in apparent magnitude, to test 
whether a possible incompleteness of the survey at bright 
and/or faint apparent magnitudes could generate a fake 
signal. To this aim we used we have 14.5 < m r < 17.5 
and we selected 370,893 objects, i.e., about 25% less than 
with less conservative constraints. We have considered 
a rectangular angular fields, with uniform coverage, in 
the SDSS internal angular coordinates (77, A) limited by 
-6° < 77 < 36° and -48° < A < 32.5°. 

To construct volume limited (VL) samples that are un- 
biased for the selection effect related to the cuts in the 
apparent magnitude, we have applied a standard proce- 
dure [33]. Firstly we compute metric distances R{z) us- 
ing the standard cosmological parameters Qm = 0.3 and 
£1a = 0.7. Secondly the galaxy absolute magnitude is de- 
termined to be M r = m r — 5-log 10 [R(z) ■ (1 + z)] — K r (z) — 
25 , where K r (z) is the K-correction. We determine the 
K r (z) term from the NYU VACG data [32]. Finally, we 
have considered two different VL samples (named VL1 
and VL2) defined by two chosen limits in absolute magni- 
tude and metric distance: for VL1 R G [100,300] Mpc/h, 
M £ [-22,-20] and for VL2 R € [200,600] Mpc/h, 
M 6 [-23, -21.5]. The number of galaxies is about 4 • 10 4 
in VL1 and 3 • 10 4 in VL2. Different cuts in absolute 
magnitude do not introduce substantial differences in the 
results presented in this paper. When more conservative 
and stringent limits in apparent magnitude are applied we 
find that in the samples with the same limits in distance, 
there are up to three times less galaxies. However the 
main results presented in this paper arc affected only in 
the fact that statistics is less robust. 

The Millennium project [38] has performed several 
cosmological simulations of standard theoretical models. 
Amount of dark matter and cosmological parameters are 
given in agreement with standard models. The dark mat- 
ter simulations have about 10 10 particles. From these 
galaxies are identified according to semi-analytics mod- 
els of galaxy formation [39]. We have cut a sample with 
exactly the same geometry as the SDSS VL1 sample and 
a sample close to the geometrical parameters of the SDSS 
VL2 applying the same absolute magnitude limits in r- 
filtcr as for the SDSS data. In the SDSS we use a redshift 
space analysis while in mock catalogs a real space one. 
The difference between the real and redshift space analy- 
sis is relevant for very small scales, i.e., r < 5 Mpc/h and 
thus does not influence results on scales of the order of 100 
Mpc/h [14,20,21]. 
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Statistical methods. — Statistical properties are de- 
termined by making averages over the whole sample vol- 
ume [2] . In doing so one implicitly assumes that a certain 
quantity measured in different regions of the sample is 
statistically stable, i.e., that fluctuations in different sub- 
regions are described by the same probability density func- 
tion (PDF). However it may happen that measurements 
in different sub-regions show systematic differences, which 
depend, for instance, on the spatial position of the spe- 
cific sub-regions. In this case the considered statistic is 
not statistically stationary in space, the PDF systemati- 
cally differs in different sub-regions and its whole-sample 
average value is not a meaningful descriptor [2] . 

In general such systematic differences may be related 
to two different possibilities: (i) that the underlying dis- 
tribution is not translationally and/or rotationally invari- 
ant; (ii) that the volumes considered are not large enough 
for fluctuations to be self-averaging [40]. On very gen- 
eral grounds, we expect the galaxy distribution to satisfy 
the condition of statistical stationarity in space to avoid 
special points or directions [1,2]. Hence the question we 
face in a finite volume analysis concerns whether it is 
large enough to obtain statistically stable results. Note 
that stationary stochastic distributions satisfy the condi- 
tion of spatial statistical isotropy and homogeneity also 
when they have zero average density in the infinite volume 
limit [2] . This condition is called the Conditional Cosmo- 
logical Principle [2] to differentiate it from the stronger 
Cosmological Principle which requires exact homogene- 
ity and deterministic rotational and translation invari- 
ancc [1,2]. 

For the case of galaxy surveys there is an intrinsic pre- 
ferred direction which is set by the radial position from the 
observer, i.e., the Earth. It is thus necessary to show that 
statistical quantities do not depend on the radial distance 
from us. To this aim, in a given sample, a simple approach 
is to determine the number N(r; R) of galaxies in spheres 
of radius r, centered on a galaxy whose distance from the 
origin is R: we call it the scale-length (SL) analysis. As we 
discuss below, this is found to be very efficient in mapping 
large scale structures which appear as large fluctuations 
of N(r; R). For instance by studying it in various angular 
slices of the SDSS samples we identify a giant filament cov- 
ering, in the largest contiguous angular area of the survey, 
more than 400 Mpc/h at R ~ 500 Mpc/h. In different sky 
directions the SL analysis reveals a variety of structures, 
showing that large density fluctuations are quite typical. 

Averaged over the whole sample the quantity N(r;R) 
gives an estimate of the average conditional number of 
galaxies in spheres of radius r. An estimator making the 
weakest a-priori-assumptions about the properties of the 
distribution outside the sample volume is [2-4] 
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Fig. 1: Left panels: From top to bottom the SL analysis for the 
SDSS sample VL1, with r = 10,20,40 Mpc/h. Right panels: 
Probability density function of N(r; R) in two non-overlapping 
sub-samples with equal volume (each half of the sample vol- 
ume) at small and large R. While for r = 10,20 Mpc/h the 
PDF (nearby sub-sample solid line and faraway sub-sample 
dotted line) is reasonably statistically stable, for r = 40 Mpc/h 
there is a clear difference. The dashed line corresponds to the 
Poisson distribution: a Gaussian function gives very good fits 
for all r. 



center-point and the number of centers M(r) varies with r 
because only those galaxies for which the sphere is fully in- 
cluded in the sample volume are considered as centers [2] . 
Even in this case, there is an intrinsic selection effect re- 
lated to the geometry of the samples, which are portions 
of spheres: when r is large only a part of the sample is 
explored by the volume average. Hence for large sphere 
radii M(r) decreases and the location of the galaxies con- 
tributing to the average in Eq.l is mostly at radial distance 
~ [Rmin + i", R m ax — r] from the radial boundaries of the 

Sample at \R m im Rmax]- 

When Eq.l scales as N(r) ~ r D and D = 3 the distri- 
bution is homogeneous, while for D < 3 it is fractal [2,24]. 

2 2 

Furthermore fluctuations 5 2 {r) = [N(r) 2 — N(r) ]/N(r) 
are small for a homogeneous distribution with any kind of 
small-amplitude correlations (S 2 (r) <C 1) and large for a 
fractal one (6 2 (r) ~ 1) [2,24]. To study fluctuations we de- 
termine the PDF of Ni(r), which is expected to converge 
to a Gaussian when r>A [2] . 

Results. — Let us now consider the VL1 sample. Here 
the SL analysis (Fig.l) detects large density fluctuations 
without a clear radial-distance dependent trend. Corre- 
spondingly the PDF has a regular shape characterized by 
a peak with a long N tail and it is sufficiently statistically 
stable in different non-overlapping sub-samples of equal 
volume. This occurs except for the largest sphere radii, 
i.e., for r >30 Mpc/h, for which the number of indepen- 
dent centers becomes too small. 
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Fig. 2: The same of Fig.l for the VL2 sample with r — 
20,40,100 Mpc/h. The fact that the signal for SDSS data 
becomes smoother when r increases, but still with a system- 
atic radial distance-trend, is due to the fact that many spheres 
overlap when r growths. While for r = 20 Mpc/h the PDF is 
reasonably statistically stable, for r = 40, 100 Mpc/h there is 
a clear difference. 

In the deeper VL2 sample we find instead a systematic 
trend of N(r;R) as a function of R (Fig. 2). Particularly 
N(r;R), for R > 300 Mpc/h, grows without any clear 
saturation, for sphere radii up to r ~ 100 Mpc/h. The 
PDF in two non-overlapping sub-samples of equal volume 
is found to differ systematically for r > 30 Mpc/h, and its 
average value moves as a function of R. This shows that 
fluctuations are not self-averaging at those scales. This be- 
havior is due to the large scale structures at scales R > 300 
Mpc/h. However at smaller scales, i.e., r < 30 Mpc/h, 
the PDF in different sub-samples is reasonably statisti- 
cally stable and similar to the one found in VL1. This 
shows that, at those scales, fluctuations are self-averaging 
because the volume average can explore different regions 
of the sample. Thus these results show that, at the largest 
scales probed, there are large density fluctuations which 
are not self-averaging because of the limited sample vol- 
ume. These determine relative variations larger than unity 
in the estimation of the average density in spheres of radius 
r = 100 Mpc/h. We thus conclude that the homogeneity 
scale must be Ao > 100 Mpc/h, the largest sphere radius 
we considered. 

Previous analyses of smaller galaxy catalogs, e.g., [4- 
6,17-19,22,33,34], considered sample averaged statistics 
without quantitatively testing whether a significant bias 
could affect the results (but see [19,22]). For instance 
the estimator of the most commonly used statistics, the 
two-point correlation function [35] , can be written as [2] 



€(r) + 1 = 



N(r, Ar) V 



(2) 



V(r,Ar) N ' 
The first ratio in the r.h.s. of Eq.2 is the average con 



ditional density, i.e., the number of galaxies in shells of 
thickness Ar averaged over the whole sample, divided by 
the volume V(r, Ar) of the shell. The second ratio in 
the r.h.s. of Eq.2 is the average density estimated in a 
sample containing N galaxies and with volume V. When 
measuring this function we implicitly assume, in a given 
sample, that: (i) fluctuations are self-averaging in differ- 
ent sub-volumes [2] (ii) the linear dimension of the sam- 
ple volume is V 1 / 3 ^> Ao [2,3], i.e., the distribution has 
reached homogeneity inside the sample volume. When 
the latter condition is not verified the £(r) analysis is bi- 
ased by systematic finite size effects even if fluctuations 
are self-averaging [2,3]. To show how non self- averaging 
fluctuations inside a given sample bias the £(r) analysis, 
we consider the estimator 



£(r; R, AR) + 1 = 



N{r, Ar) 



V(r*) 



V(r,Ar) N(r*;R,AR) 



(3) 



where the second ratio on the r.h.s. is the density of points 
in spheres of radius r* averaged over the centers lying in a 
shell of thickness AR around the radial distance R. If the 
distribution is homogeneous, i.e., r* > Ao, and statistically 
stationary, Eq.3 should be statistically independent on the 
range of radial distances (R, AR) considered. For instance 
we consider, in the VL2 sample, AR = 40 Mpc/h and R = 
240 Mpc/h or R = 520 Mpc/h, with r* > 50 Mpc/h. We 
thus find large variations in the amplitude of £(r) (Fig. 3). 
This is simply an artifact generated by the large density 
fluctuations on scales of the order of the sample sizes. The 
results that the estimator Eq.2 (or others based on pair- 
counting [2,36]) has nearly the same amplitude in different 
samples, e.g., [7-9,33,34], despite the large fluctuations of 
N(r;R), are simply explained by the fact that £(r) is a 
ratio between the local conditional density and the sample 
average density: both vary in the same way when the 
radial distance is changed and thus the amplitude is nearly 
constant. 

On the other hand Eq.l, averaged over the volumes 
where the PDF has a statistically stable shape, shows in 
both considered samples a power-law behavior for r < 30 
Mpc/h corresponding to a fractal dimension D — 2.1 ±0.1 
in agreement with [4,6,17,19,22] (Fig. 3). Due to the non 
self-averaging nature of fluctuations at larger scales, i.e., 
due to limited volumes, we are not able to determine cor- 
relations for r > 30 Mpc/h. 

Discussion. — According to standard models of cos- 
mological structure formation, gravitational clustering 
gives rise to non-linear perturbations from homogeneous 
initial conditions in the early universe [1]. If the initial 
amplitude of fluctuations is normalized to the anisotropies 
of the Cosmic Microwave Background Radiation (CMBR) 
[37], then the homogeneity scale is about A™ = 10 Mpc/h 
[23], i.e., twice the value at which £(r) = 1 [1,2]. 

Indeed in mock galaxy catalogs generated from N- 
body simulations of standard cosmological models [38,39], 
N(r;R) does not show, for r > A™, large fluctuations 
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Fig. 3: Upper Panel: The sample average conditional num- 
ber of galaxies (Eq.l) for the SDSS VL1 and VL2 samples 
up to r < 30 Mpc/h. The best fit slope (soldi line) gives 
D = 2.1 ± 0.1. The difference in amplitude between the two 
samples is simply ascribed to a luminosity selection effect [2], 
as VL1 contains fainter galaxies than VL2. The amplitude of 
the mock VL1 and VL2 samples has been rescaled by the same 
arbitrary factor for seek of clarity. Bottom Panel: Standard 
two-point correlation function in the VL2 sample estimated by 
Eq.3: the sample average density is computed in spheres of 
radius r* and considering all center-points lying in a bin of 
thickness AR = 40 Mpc/h centered at different radial distance 
R. The case r* = 80 Mpc/h and R = 520 Mpc/h gives the 
average over the whole sample, i.e., Eq.2, and it coincides with 
the estimation of [33] in a similar sample. The flat tail of f (r), 
for r* = 40,80 Mpc/h and R = 300 Mpc/h, reflects the inho- 
mogeneities with strong correlations at large scales. The stars 
and triangles correspond to the behavior of the mock VL1 and 
VL2 samples: in this case the amplitude is statistically stable 
and thus meaningful. The small amplitude difference in this 
case is ascribed to the different selection in luminosity [38]. 

or systomatic trends as a function of R (Fig. 4). Because 
in these artificial catalogs fluctuations are small and self- 
averaging, whole-sample averaged statistics are meaning- 
ful at all scales. From the N(r) analysis we find, differ- 
ently from the real galaxy data, that N(r) ~ r D with 
D = 1.1 ± 0.1 for r < Xfi 1 and D = 3 for r > A™ (Fig.3). 
Correspondingly the PDF rapidly converges to a Gaussian 
for r > A™. The £(r) function, estimated by Eq.2 or Eq.3, 
has a statistically stable amplitude. The residual small 
amplitude difference between the mock VL1 and VL2 sam- 
ples is ascribed to a different selection in luminosity [38]. 
The amplitude of £(r) is equal to unity at r ps 6 Mpc/h in 
agreement with [38,39]. In addition its shape reasonably 
agrees with the standard estimation of £(r) from galaxy 
catalogs [7-9,33,34]. However, as discussed above, the lat- 
ter estimation is biased by systematic effects making the 
agreement fortuitous. 

Standard theoretical models predict that, for r > A™, 
the small fluctuations in the early universe are linearly am- 
plified by gravitational clustering. Therefore for r > X™ 
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Fig. 4: The same of fig. 1 for mock sample VL2, with r = 
20,40,100 Mpc/h. Although for r = 20 Mpc/h fluctuations 
are still important, they rapidly become small for larger radii. 
The PDF is statistically stable for all r and R. The solid line 
corresponds to a Gaussian fit. 

the shape of the theoretical £ m (r) must be the same as 
the initial one [1] . This is characterized by a length-scale 
r c , where £, m ( r c) — 0, which is fixed by the physics of 
the early universe and estimated from CMBR anisotropics 
to be r c w 100 Mpc/h [36,37]. For r > r c , f™(r) be- 
comes negative, corresponding to super-homogeneous cor- 
relations characterized by the most rapid possible decay 
of fluctuations [41,42]. This theoretical framework ap- 
plies to the whole mass distribution, where dark matter is 
supposed to provide with the main contribution. Galax- 
ies would form on the largest peaks of the density field. 
Standard models of galaxy formation describe this physi- 
cal phenomenon as a selection mechanism [44] . This leaves 
unperturbed the scale r c and slightly changes A™ [36,43]. 
We find A > 100 Mpc/h > r c > Aft 1 . This raises a 
fundamental inconsistency for the relation between galaxy 
structures and CMBR anisotropics as no physical mecha- 
nism is known, which by sampling a super-homogeneous 
density field transforms it into a strongly inhomogeneous 
one [2,36,43]. 

Conclusion. — In summary, by applying the SL anal- 
ysis to the newest SDSS galaxy samples, we measure large 
density fluctuations of spatial extension limited by sam- 
ple sizes. At scales r < 30 Mpc/h we detect statistically 
stable fractal correlations with D = 2.1 ± 0.1. On larger 
scales, r > 30 Mpc/h, we find that the galaxy distribu- 
tion is strongly inhomogeneous and fluctuations are not 
self-averaging in the samples considered. This situation is 
compatible with fractal power-law correlations extending 
to such length-scales but incompatible with homogene- 
ity at A < 100 Mpc/h. Indeed, in a portion of a frac- 
tal, large structures are expected to be present at any 
scale, fluctuations being self-averaging only if the sam- 
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pie volume is large enough [2]. These results have impor- 
tant consequences on the theoretical interpretation of the 
large scale universe, where models, normalized to CMBR 
anisotropies, predict there is not enough time to form 
structures with relative density fluctuations larger than 
unity on scales larger than A™ w lOMpc/h [1,38]. This 
length scale is more than ten times smaller than our lower 
limit to Ao. Indeed the latter is of the order of the scale 
r c where theoretical model predict matter distribution to 
have negative correlations, a situation which is in contrast 
with the results from the data analyzed here. Thus the 
large scale inhomogeneities detected in the SDSS samples 
are incompatible with the predictions of standard theo- 
retical models relating the early universe physics, with 
CMBR normalization, to structures in the present uni- 
verse. Moreover we found that for r < A™, mock galaxy 
catalogs have different correlations from real galaxy data, 
i.e., D = l.liO.l instead of D = 2.1±0.1. Thus structures 
generated by N-body simulations arc intrinsically different 
from observed ones. 

Recently three dimensional maps of dark matter dis- 
tribution from weak lensing observations have been pub- 
lished [45]. Dark matter is observed to trace the same 
structures as galaxies. Thus the whole matter distribu- 
tion would be inhomogeneous on scales larger than 100 
Mpc/h. This has a great impact on the whole theoretical 
framework for the physical understanding of the large scale 
universe. For instance it may imply a new type of evolu- 
tion scenario within an open Fricdmann model [46] or new 
types of spatial averaging of the Einstein equations [47,48] 
which relate the observed inhomogeneities to the apparent 
acceleration measured from supernovae observations [49] . 

The determination of statistical properties of those very 
large structures, which we detected but could not analyze 
in detail, should be possible when sample volumes become 
large enough so that the corresponding fluctuations will 
be self-averaging on scales larger than the ones studied 
here. The application of the SL analysis to the forthcom- 
ing galaxy redshift surveys, like the complete SDSS [15], 
thus represents an important task. 
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